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To  Compute  the  conjugate  depth  at  a  hydrau 
triangular,  or  rectangular  open  channel^ 

References:  (l)  Chow,  Ven  Te,  Open  Channel  Hydraulics,  McGraw-Hill,  1959 
PP  53-56 

(2)  Schied,  Francis,  Numerical  Analysis,  Schaum's  Outline  Series, 
McGraw-Hill,  1968,  pp  315-317* 


I.  PHOQRAM  SPCCIFtCATtOMt 

SEE  FOLLOWING  PAGE 


The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is  part  of 
a  Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS)*  The 
program  consists  of  a  main  program  and  three  subroutines.  The  main  program 
handles  all  I/O  requirements  with  the  subroutines  handling  the  computation.^ 


».  IQOIPMCMT  DCTAIlT 


rhe  program  was  developed  and  is  operational  on  the  WES  G635>  Vicksburg,  MS. 
It  is  also  operational  on  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 


«.  IMPUT.WTWr - 

The  required  inputs  which  are  entered  at  execute  time  are:  discharge,  cfs, 
bottom  width,  ft,  (0  if  tria  sect.),  side  slope  (0  if  rect  sect),  number 
of  flow  depths  and  the  flow  depths,  ft.  Output  includes  the  type  of  cross 
section  and  the  given  data,  plus  the  conjugate  depth,  ft. 


ADDITIONAL  NCNANK* 


[-Complete  documentation  of  this  program  is  available  from  the 
Engineer  Computer  Programs  Library,  Technical  Information  Center,  WES. 
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B.  PROGRAM  SPECIFICATIONS: 

Language :  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements:  The  run  command 

RUN  WESLIB/C0RPS/H6123,  R 

and  the  inputs  defined  in  (E). 

Method  of  Analysis:  Solves  the  momentum  equation  for  the  conjugate 
depth  hy  application  of  Newton* s  Tangent  Method. 

Core  Requirements  G635 :  11  K  words 

External  Storage:  None 

Restrictions :  Velocity  distribution  is  uniform  across  the  flow 
section  and  the  invert  slope  is  small  (i.e.,  <  10  degrees). 

General  Equation: 


2  2 

Vi  ♦  Iq’ =  Vs  *  isj 

where  Q  is  discharge  (cfs),  g  is  acceleration  of  gravity 
2 

(32.2  ft/sec  ),  A^  and  A^  are  the  respective  cross-section  areas 
(ft  )  for  flow  depth  and  conjugate  depth,  and  y^  and  y^  are  the 
distances  (ft)  of  the  centroids  of  the  respective  areas  A^  and  A^ 
below  the  water  surface. 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer 
Program  Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-R0-6AI 

2.  TITLE:  H6123  -  Conjugate  Depth  in  Either  a  Trapezoidal,  Triangular, 
or  Rectangular  Open  Channel. 

3.  REVISION  LOG:  N/A 

4.  PURPOSE  OF  PROGRAM:  To  compute  the  conjugate  depth  at  a  hydraulic 
jump  in  either  a  trapezoidal,  triangular,  or  rectangular  open  channel. 
References: 

a.  Chow,  Ven  Te,  Open  Channel  Hydraulics,  McGraw-Hill,  1959* 

pp  53-56. 

b.  Schied,  Francis,  Numerical  Analysis,  Schauta1  s  Outline  Series, 
McGraw-Hill,  1968,  pp  315-317- 

5-  STEP  SOLUTION: 

a.  Enter  inputs: 

(1)  Discharge  (Q),  cfs 

(2)  Channel  bottom  width  (b),  ft  (0  if  triangular  section) 

(3)  Channel  side  slope  (z),  expressed  as  the  cotangent  (0  if 
rectangular  section) 

(4)  Flow  depth  (y^), 

b.  Computational  steps: 

(1)  Critical  depth  (yQ)  using  subroutine  H6l4o 

(2)  Specific  force  (F^)  at  y^  using  subroutine  H6124 

(3)  Determine  trial  conjugate  depth  (yt) 


It  yc  <  y±t  then  yt  -  j- 


SPECIFIC  FORCE  CURVE 


If  yc  =  yx»  then  the  conjugate  depth  (yg)  is  equal  to 
yc  ;  yg  is  returned  to  the  main  program.  If  yfi  >  y^  » 


SPECIFIC  FORCE  CURVE 


(U)  Specific  force  (F  )  and  area  (A  )  at  y.  using  sub¬ 
routine  H612U 

(5)  Define  function  F(y.)  to  be  used  for  Newton's  Tangent 
Method 


F(yt)  =  Ft  -  Fx 

(6)  First  derivative  of  F(yt)  with  respect  to  yt 


.  Q2(b  +  2zy  ) 

F(yJ  m  K  - 3 — 
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(7)  Conjugate  depth  (yp)  is  calculated  using  Newton's  Tangent 
Method 

y 

(8)  If  y  £  0.0  ,  then  y.  =  —  ,  where  m  *  2,3,  • . . , depending 

tL  w  2m 

on  the  number  of  times  the  condition  y2  £  0.0  occurs. 
Transfer  is  to  step  (U)  and  the  process  is  repeated  until 
yg  >  0.0.  Step  (9)  is  then  entered. 

(9)  If  |yg  -  yt|  >  .001  ,  then  y^  =  y^  and  the  procedure 

returns  to  step  (U)  and  is  repeated  until  | y g  —  y^ |  <  .001. 
The  value  of  i-s  returned  to  the  main  program. 

c.  The  given  data,  plus  the  conjugate  depth,  are  printed. 

6.  ACCURACY :  Governed  by  accuracy  of  input  data;  conjugate  depth  is 
computed  to  +_  .001  ft. 

7.  REMARKS:  Velocity  distribution  is  uniform  across  flow  section  and 
the  invert  slope  is  small  (i.e.,  <  10  degrees). 
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3.  EQUIPMENT  AND  OPERATING  SYSTEM:  The  program  was  developed  on  a 
G635  time-share  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype.  It  is  now  operational  on  the  WES  G635, 
Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

1*.  INPUT  REQUIREMENTS:  The  required  inputs  are  entered  via  the  user's 
time-share  terminal  device  in  free  field  format.  All  input  cues  and 
reads  are  performed  in  the  main  program.  The  subroutines  handle  the 
computations.  Since  computations  are  done  in  the  subroutines,  and  sub¬ 
routines  H6lU0  and  H6l24  are  subordinate  to  H6l23»  the  necessary  inputs 
to  subroutine  H6123  ere  passed  via  the  CALL  statement.  The  calling 
sequence  is: 

CALL  H6123  ( ar , . . . ,arg^ ) 

where : 

arg^  -  discharge,  cfs 
argg  -  bottom  width,  ft 
arg^  -  side  slope  ( cotangent ) 
arg^  -  flow  depth,  ft 
arg^  -  conjugate  depth,  ft 

All  arguments  are  floating  point.  Arguments  1-1*  are  the  inputs  and 
argument  5  is  the  output. 

5 .  SECONDARY  STORAGE  INPUT:  None 

6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  H6123. 
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DISCH 

-  discharge,  cfs 

(arg^ 

BWID 

-  bottom  width,  ft 

(argg) 

SSLOPE 

-  side  slope  (cotangent) 

(arg^) 

N 

-  number  of  flow  depths,  integer 

Y 

-  flow  depth,  ft 

(arg^) 

7- 

OUTPUT 

DATA  DESCRIPTION:  The  following  name  is  used  for  the  output 

variable  in  program  H6123 

DEPTH2 

-  conjugate  depth,  ft 

(arg^) 

8. 

PROGRAM  ERROR  MESSAGES:  None 

9. 

VARIABLE  DEFINITIONS: 

AREA 

p 

-  area  at  trial  conjugate  depth,  ft 

BWID 

-  bottom  width,  ft  (0  if  tria  sect) 

CRIT 

-  critical  depth,  ft;  returned  from  subroutine  H6lU0 

N 

-  number  of  flow  depths,  0  <  N  <  25 

Y 

-  flow  depth,  ft;  dimensioned  max  25 

DEPTH1 

-  flow  depth,  ft;  equal  to  Y(l)  for  I  *  1,...,N 

DEPTH2 

-  conjugate  depth,  ft 

DEPTHT 

-  trial  conjugate  depth,  ft 

DISCH 

-  discharge,  cfs 

F 

-  function  of  trial  conjugate  depth  used  for  Newton's 

Tangent  Method,  ft^;  the  value  of  the  function  is  a 
specific  force  difference  between  that  at  trial  conjugate 
depth  and  flow  depth 

FP 

-  first  derivative  of  function  (F)  with  respect  to  trial 
conjugate  depth,  ft2 
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3 

FORCE  -  specific  force,  ft  ;  returned  from  subroutine  H6l24 

FI  -  working  storage  to  hold  FORCE  at  a  constant  when  FORCE  is 

returned  from  H6124  as  the  specific  force  at  flow  depth, 
ft3 

p 

G  -  acceleration  of  gravity,  32.2  ft/sec 

HFILE  -  five  character  name  of  program;  passed  to  WESLIB  count 
routine  HACCT 

LQZ  -  equal  1,  execute  all  input  cues  and  reads;  equal  2,  call 
WESLIB  routine  RERUN  and  enter  only  desired  inputs 

LQX  -  equal  1,  print  instructions  from  RERUN;  equal  3,  no  print 

JKL  -  direct  return  from  RERUN  to  desired  input  read 

KKK  -  total  number  of  inputs;  passed  to  RERUN 

M  -  power  of  2  to  converge  Newton's  Tangent  Method  if  conju¬ 

gate  depth  0 

TYPE  -  character;  type  of  flow  cross  section;  trap,  tria,  or  rect 

ZZZZZ  -  character;  equal  RE,  rerun;  equal  ST,  stop 
10.  EXAMPLE  CASE:  Compute  the  conjugate  depth  for  6  given  flow  depths 
in  a  rectangular  channel . 


Flow  depths  (Y(I),  for  I  *  1,  N)  »  1.00,  1.43,  2.50,  4.00, 
5.27,  and  6.00  ft 
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b.  Output : 


INPUT  H6 123-CONJUGATE  DEPTH-TRAP, RECT, TRIA-OPEN  CHANNEL 

AA-ENTER  DISCHARGE# CFS . 

=  100 

AB-ENTER  CHAN  BOTTOM  WIDTH, FT  (0  IF  TRIA  SECT) . 

=  3.51 

AC-ENTER  CHAN  SIDE  SLOPE,  (0  IF  RECT  SECT). 

=  0 

AD-ENTER  THE  NUMBER  OF  DEPTH(S)  FOR  WIHICH  THE  CONJUGATE  DEPTH  IS  TO  BE 
CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS. 

=  6 

AE-ENTER  THE  6  DEPTH(S)  SEPARATED  BY  COMMAS. 

=1,1.43,2.5,4,5.27,6.5 

OUTPUT  H6123-C0NJUGATE  DEPTH(S)  FOR  A  RECT  OPEN  CHANNEL 

DISCHARGE  =  100.00  CFS 

BOTTOM  WIDTH  =  3.51  FT 


SIDE  SLOPE 

— 

FLOW 

CONJUGA 

DEPTH 

DEPTH 

(FT) 

(FT) 

1.00 

6.62 

1.43 

5.27 

2.50 

3.41 

4.00 

2.07 

5.27 

1.43 

6.50 

1.03 

ENTER  RERUN  OR  STOP 
=STOP 
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REP:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H6123  -  Conjugate  Depth  in  Either  a  Trapezoidal,  Triangular, 
or  Rectangular  Open  Channel. 

3.  PROGRAM  SOURCE  LISTINGS:  See  pages  12-lU 

4.  NUMERICAL  AND  LOGICAL  ANALYSIS:  Solves  the  momentum  equation  for 
the  conjugate  depth  by  application  of  Newton's  Tangent  Method. 

5-  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  The  following  subroutines 
are  used  in  program  H6123. 

a.  SUBROUTINE  H6l24(DISCH,  BWID,  SSLOPE,  DEPTH,  AREA,  FORCE) 

The  subroutine  statement  for  H6124  as  documented  in  program 
H6124  is: 

SUBROUTINE  H6l24(DISCH,  BWID,  SSLOPE,  DEPTH,  ENERGY,  FORCE) 

The  argument  list  was  changed  to  (a)  to  facilitate  the  solution 
of  program  H6123.  This  change  has  no  effect  on  H6124  as 
documented. 

b.  SUBROUTINE  H6l4o(DISCH,  WIDTH,  SSLOPE,  DEPTH,  CRIVEL) 

Complete  documentation  of  these  subroutines  is  available  from 
the  Engineer  Computer  Program  Library,  Technical  Information 
Center,  WES. 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational,  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
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The  source  listing  on  page  12  contains  the  first  line  run  command  and 
brief  for  H6123.  This  first  line  run  command  runs  the  binary  H6123B  of 
the  source  listing  on  pages  13-lU  (Fortran  source  of  H6123)  and  attaches 
the  WESLIB  routines  HACCT  and  RERUN. 
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OOOlKtRUN  WESLIB/CORPS/H6123B.R.WESL IB/RERUN, RjUESLIB/HACCT.R 

0800  63THJ5  PROGRAM  COMPUTES  CONJUGATE  DEPTHS  IN  EITHER  A  TRAPEZOIDAL. 
0810  62TRIANGULAR,  OR  RECTANGULAR  CHANNEL.  INPUTS  REQUIRED  ARE  CHAN- 
0820  59NEL  BOTTOM  UIDTH  IN  FTtZERO  IF  TRIANGULAR  SECTION).  CHANNEL 
0830  61SIDE  SLOPE  EXPRESSED  AS  THE  COTANGENT  OF  THE  ACUTE  ANGLE  WITH 

0840  63RESPECT  TO  THE  H0RIZ0NTAV<^|R?».jriiSliTTlJpUrT5FHEnlTAN>TYPEDOFHE 
0850  61FL0M  DEPTH(S)  IN  FT.  OUTPUT  INCLUDES  THE  GIVEN  DATA,  TYPE  OF 

0860  42CR0SS-SECTI0N,  AND  THE  CONJUGATE  DEPTH(S) . 


0999X06FINISH 


12 


H6123 


00001*»RUN  *=;/C0RPS/H6123B(N0G0) 
lOOOO  CHARACTER  TYPE»4,HFILE*5 
10010  DIMENSION  DEPTH(25) 

10020  HFILE=5HH6123 

10030  LQZ=1;LQX=1 

10040  13000  CALL  HACCT(HFIIE) 

10050  PRINT  15111 

10060  15111  FORMAT (/"INPUT  H6 123-CONJUGATE  DEPTH-TRAP, RECT.TRIA-OPEN 
100650  CHANNEL"//) 

10070  GO  T0(15003, 15023), LQZ 
10080  15003  PRINT  15004 

10090  15004  FORMAT ("A A- ENTER  DISCHARGE, CFS.") 

10100  15005  READ, DISCH 
10110  GO  TOU5007, 15023), LQZ 
10120  15007  PRINT  15008 

10130  15008  FORMAT ( "AB-ENTER  CHAN  BOTTOM  WIDTH, FT  (0  IF  TRIA  SECT).") 

10140  15009  READ, BUID 

10150  GO  TOClSOll, 15023), LQZ 

10160  15011  PRINT  15012 

10170  15012  FORMAT ( "AC-ENTER  CHAN  SIDE  SLOPE,  (0  IF  RECT  SECT).") 

10180  15013  READ, SSLOPE 
10190  GO  T0(15015, 15023), LQZ 
10200  15015  PRINT  15016 

10210  15016  FORMAT("AD-ENTER  THE  NUMBER  OF  DEPTH(S)  FOR  UIHICH  THE  CONJUG 
102150ATE  DEPTH  IS  TO  BEV’CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS.") 

10220  15017  READ,N 

10230  GO  TO(15019, 15023), LQZ 

10240  15019  PRINT  15020, N 

10250  15020  FORMAT ("AE- ENTER  THE  ",I2."  DEPTH(S)  SEPARATED  BY  COMMAS.") 
10260  15021  READ, (DEPTH(I), 1=1, N) 

10270  GO  TOU5026, 15023), LQZ 

10280  15023  KKK=5 

10290  CALL  RERUN (KKK,LQX,JKL) 

10300  GO  TOU5005, 15009, 15013, 15017, 15021, 15026), JKL 

10310  15026  TYPE=4HTRAP 

10320  IF(BUID)  15029,15028,15029 

10330  15028  TYPE=4HTRIA 

10340  15029  IF(SSLOPE)  15031,15030,15031 

10350  15030  TYPE=4HRECT 

10360  15031  PRINT  15032, TYPE. DISCH, BUID, SSLOPE 

10370  15032  FORMAT (/"OUTPUT  H6 123-CONJUGATE  DEPTH(S)  FOR  A  ",A4,"  OPEN 
103800  CHANNEL"// "DISCHARGE  =  ",F12.2,"  CFS"/"BOTTOM  WIDTH  =  ",F9.2,"  FT 
103900"/"SIDE  SLOPE  =  ", Fll .2, "H* 1V"//4X, "FLOW", 5X, "C0NJUGATE"/4X, "DEPTH 
104000", 6X,"DEPTH"/4X,"( FT)", 7X,"( FT)") 

10410  DO  15039  1=1, N 
10420  DEPTII1=DEPTH(I) 

10430  CALL  H6123(DISCH,BWID. SSLOPE, DEPTH1 ,DEPTH2) 

10440  15039  PRINT  15040, DEPTH1 .DEPTH2 
10450  15040  F0RMAT(F8.2,F11 .2) 

10440  PRINT,"  " 
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10470  LQZ=2 

10480  CHARACTER  ZZZZZZ#2 

10490  16000  PRINT ,  "ENTER  RERUN  OR  STOP" 

10500  READ  16001 »  ZZZZZZ 
10510  16001  F0RMAKA2) 

10520  IF (ZZZZZZ . EQ . 2MRE )  GO  TO  15000 
.10530  IF (ZZZZZZ . EQ . 2HST)  GO  TO  20000 
10540  PRINT, "ERROR  ###  RETYPE" 

10550  GO  TO  16000 
10560  20000  STOP ; END 

20000  SUBROUTINE  H6123< DISCH, BWID, SSLOPE, DEPTH1 ,DEPTH2) 

20010  G=32 . 2 ; M=2 

20020  CALL  H6 1 40 (DISCH , BWID, SSL0PE, CRIT , CV ) 

20030  CALL  H6 124 (DISCH, BWID, SSLOPE, DEPTH 1 , AREA, FORCE) 

20040  Fl=FORCE 

20050  IF(CRIT-DEPTHl)  20060,20070,20080 
20060  20060  DEPTHT=CRIT/N ; N=2*N ; GO  TO  20090 
20070  20070  DEPTH2=CRJ T ; RETURN 
20080  20080  DEPTHT  =CRIT##2/DEPTH1 

20  0  90  20090  CALL  H6 1 24 ( DISCH, BU ID, SSL OPE, DEPTHT, AREA, FORCE) 

20100  F=F0RCE-F1 

20110  FP=AREA-DISCH##2#(BWID+2.#SSLOPE#DEPTHT)/(G*AREA*#2) 

20120  DEPTH2=DEPTHT-F/FP 

20130  IF( DEPTH2 . LE. 0 . 0 )  GO  TO  20060 

20140  IF(ABS(DEPTH2~DEPTHT).LE. .001)  RETURN 

20150  DEPTHT=DEPTH2;G0  TO  20090 

20160  END 

30000  SUBROUTINE  H6124 ( DISCH , BWID, SSLOPE, DEPTH, AREA , FORCE) 

30010  G=32.2 

30020  AREA - DEPTH# ( BWID*SSLOPE#DEPTH) 

30030  EHERGY=DEPTH+DISCH##2/( 2 . #G#AREA##2) 

30040  FORCE=BWID# DEPTH# #2/2 . + SSL OP E# DEPTH## 3/3. +DISCH# #2/(0# AREA) 
30050  RETURN 
30060  END 

40000  SUBROUTINE  H6140 ( DISCH , WIDTH , SSLOPE, DEPTH, CRIVEL ) 

40010  G-32 . 2 

40020  IF(UIDTH)  30,10,30 

40030  30  DEPTH=( (DI$CH##2)/(G#(WIDTH##2) ) )##(1 ./3 .  ) 

40040  GO  TO  20 

40050  10  DEPTH=((2.#DI$CH##2)/(G#$SlOPE##2))##(l./5.) 

40060  GO  TO  100 

40070  20  IF( SSLOPE)  2 5,100,25 

40080  25  A=DISCH##2#(WIDTH+2.#SSL0PE#DEPTH) 

40090  B=G#DEPTH##3#CWIDTH>SSLOPE#DEPTH)##3 
40100  DIFF=B-A 

40110  IF(DIFF.LT. 0.0001)  GO  TO  100 
40120  DEPTH  =  DEPTH-0.01 
40130  GO  TO  25 
40140  100  CONTINUE 

40150  AREA=(UIDTH+SSLQPE#DEPTH)#DEPTH 
40160  CRIVEL  -  DISCH/AREA 
40170  RETURN 
40180  END 
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